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Abstract: With the newly available data sets on hadronic r decays from the 5-factories 
BABAR and BELLE, and future data from BESIII, precise information on the decay 
distributions will soon become available. This calls for an improvement of the decay spectra 
also on the theoretical side. In this work, the distribution function for the decay r — ► v T K-n 
will be presented with the relevant Kir vector and scalar form factors being calculated in the 
framework of the resonance chiral theory, also taking into account additional constraints 
from dispersion relations and short-distances. As a by-product the slope and curvature of 
the vector form factor F±(s) are determined to be A+ = 25.6 • 1(T 3 and A+ = 1.31 • 1(T 3 
respectively. From our approach it appears that it should be possible to obtain information 
on the low lying scalar A"q(800) as well as the second vector If* (1410) resonances from the 
r decay data. In particular, the exclusive branching fraction of the scalar component is 
found to be B[r -> ^(^)s- wav c] = (3.88 ± 0.19) ■ 10" 4 . 
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1 Introduction 



Already more than a decade ago it was realised that the hadronic decays of the r lepton 
could serve as an ideal system to study low-energy QCD under rather clean conditions 
[1-5]. In the following years, detailed investigations of the r hadronic width as well as 
invariant mass distributions have allowed to determine many QCD parameters, a most 
prominent example being the QCD coupling a s . Especially the experimental separation 
of the Cabibbo-allowed decays and Cabibbo-suppressed modes into strange particles [6-8] 
opened a means to also determine the quark- mixing matrix element |V^ S | [9-11] as well 
as the mass of the strange quark [12-19], additional fundamental parameters within the 
Standard Model, from the r strange spectral function. 

The dominant contribution to the Cabibbo-suppressed r decay rate is due to the decay 
t — > u T K7r. The corresponding distribution function has been measured experimentally in 
the past by ALEPH [8] and OPAL [7] . With the large data sets on hadronic r decays from 
the B-factories BABAR and BELLE, which are currently under investigation, and good 
prospects for additional data from BESIII in the future, a refined theoretical understanding 
of the spectral functions is called for. For the decay in question, the general expression for 
the differential decay distribution takes the form [20] 

dT K7T G 2 F \V US \ 2 M? / s 
dy/s 32tt 3 s V M? 

where we have assumed isospin invariance and have summed over the two possible decays 
t~ — > v T K 7r~ and r~ — > v T K~ir°, with the individual decay channels contributing in the 
ratio 2 : 1 respectively. In this expression, F+*(s) and F^ 7r (s) are the vector and scalar Kn 
form factors respectively which will be explicated in more detail below, A KlT = M\ — M 2 , 
and qKn is the Kaon momentum in the rest frame of the hadronic system, 

q Kn (s) = ^j=^J (s - (M K + M n )^ (s - (M K - M^j ■ e(s - (M K + ikQ 2 ) . (2) 

By far the dominant contribution to the decay distribution originates from the if* (892) 
meson. In the next section, an effective description of this contribution to the vector form 
factor F^ n (s) will be presented in the framework of chiral perturbation theory with reso- 
nances (RxPT) [21,22], quite analogously to a similar description of the pion form factor 
given in refs. [23-25]. To also include a second vector resonance, the effective chiral descrip- 
tion can be straightforwardly augmented by an additional if* (1410) meson. Finally, the 
scalar Kir form factor F i ^ 7r (s) has very recently been updated in ref. [26]. In section 3, we 
shall present our main results for the distribution function dT Kn / dy/s and total decay rates 
and investigate the influence of the different vector and scalar form factor contributions. 
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2 The vector form factor F^ n (s) 

A theoretical representation of the vector form factor F+ n (s), which is based on funda- 
mental principles, can be provided by an effective field theory description, in complete 
analogy to the description of the pion form factor presented in refs. [23-25] l . This ap- 
proach employs our present knowledge on effective hadronic theories, short-distance QCD, 
the large- Nc expansion as well as analyticity and unitarity. For the pion form factor the 
resulting expressions provide a very good description of the experimental data [23-25]. 

Close to s equal zero, F+ n (s) is well described by the %PT result at one loop [29], 
which takes the form: 



H K „(s) + H K ,(s) . (3) 



The one-loop function H(s) is related to the corresponding function H(s) of [29] by H(s) = 
H(s) - 2L r 9 s/(3F%) = [sM r (s) - L(s)]/F%. Explicit expressions for M r (s) and L(s) can 
be found in ref. [30]. The vector form factor is an analytic function in the complex s 
plane, except for a cut along the positive real axis, starting at s Kn = (M K + M w ) 2 , where 
its imaginary part develops a discontinuity. In the elastic region below roughly 1 GeV, 
F+*(s) admits the well-known Omnes representation [31] 



F««( 8 ) = P( S )exp^ Jds'^ 



S 1/2 (s') 

6l [S) 1 (4) 



s - iO) J ' 



with P(s) being a real polynomial to take care of the zeros of F^ 7T (s) for finite s and <5^ 2 (s) 
is the P-wave 1 = 1/2 elastic Kit phase shift. 

Precisely following the approach of ref. [23], and matching the Omnes formula (4) with 
the %PT calculation of F^ 7r (s) in the presence of vector resonances [21], one finds the 
following representation of the form factor F^(s): 



te ^Re[H K7T (s)+H Kv (s)] 

M 2 K * - s - iM K *T K *{s) 



F K7r (s) = K - - (5) 



The one- loop function H(s) depends on the chiral scale /x, and in eq. (5), this scale should be 
taken as fi = M^*. In ref. [32], the off-shell width of a vector resonance was defined through 
the two-point vector current correlator, performing a Dyson- Schwinger resummation within 
R%PT [21,22]. Following this scheme the energy-dependent width r^*(s) is found to be 

r ^* (s) = ^I^K (s)+ ^ ( 4 (6) 



1 For an alternative dispersive approach to the pion form factor see also refs. [27,28]. 
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where Gy is the chiral vector coupling which appears in the framework of with 
resonances [21], the phase space function cr K7T (s) is given by (Jk-k{s) = 2q K7T (s)/ \/s, and 
&Kr)(s) follows analogously with the replacement M n — > M v . Re-expanding eq. (5) in s 
and comparing with eq. (3), one reproduces the short-distance constraint for the vector 
coupling G v = F n /V2 [22] which guarantees a vanishing form factor at s to infinity, as 
well as the lowest-resonance estimate 2 for the 0{p A ) chiral coupling [21] 

L r 9 (M K *) = JjL = 5.34 • 10- 3 , (7) 

where F w = 92.4 MeV and the average mass of the charged and neutral if* (892), M K * = 
893.9 MeV have been used. This result is in very good agreement to a recent determination 
of L\ from the pion form factor [34] which found L r 9 (M K *) = (5.69 ± 0.41) • 1(T 3 . 
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Figure 1: Ktc scattering P-wave phase-shift data of ref. [35], together with our fit described 
in the text. 

1 /2 

From eq. (5), one obtains a description of the P-wave Kn phase shift 5/ (s): 

' Mk*Tk*(s) 



Sy 2 (s) = arctan^-' 



Improving the if*-meson width T K *(s) with the pertinent Blatt-Weisskopf barrier factor 
Di(rq Kw (Mj c „)) / D 1 (rq K7T (s)) [36], where Di(x) = 1 + x 2 and r is the interaction radius, we 

2 A recent one-loop calculation of the pion form factor in RxPT [33] has found that the NLO corrections 
to this lowest-resonance approximation are very small. 
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can perform a fit of eq. (8) to the phase-shift data of Aston et al. [35]. The fit parameters 
M K *, Gy and r are found to be 

M K * = 895.9 ± 0.5 MeV, Gy = 64.6 ± 0.4 MeV, r = 3.5 ± 0.6 GeV~ x , (9) 

where the given errors are purely statistical and the corresponding fit is displayed in figure 1. 

We observe that the obtained K* mass is in complete agreement with the fit re- 
sult of ref. [35], and also with the PDG average for the mass of the neutral K* meson 
Mk* = 896.1 ± 0.3 MeV [37], the relevant channel in this case. Our result for Gy is nicely 
compatible with the short-distance constraint Gy = F n /y/2 ~ 65.3 MeV found in ref. [22]. 
Furthermore, our fit result for r and our % 2 /NDF m 18.5/12 are similar to the ones ob- 
tained by Aston et al. [35]. However, the result for Gy found above would correspond to 
a K* width 

T K * = T K *{M 2 K *) = 55.8 ±0.8 MeV, (10) 

which is about 10% larger than the finding of [35] and the PDG average Tk* = 50.7 ± 
0.6 MeV [37]. At the origin of this discrepancy may lie the role played by a heavier vector 
resonance, namely the if* (1410), on which we shall comment further below. As concluded 
in ref. [25], Gy goes down a few percent when a second multiplet is included into the 
analysis of the elastic tttt scattering amplitude, providing Gy = 61.9 ± 1.5 MeV [25] which 
would result in Tk* = 51.2 ± 0.7 MeV, in excellent agreement with the PDG average. 

Next, we can compare the low-energy expansion of the vector form factor (5) with recent 
experimental measurements of its slope and curvature. Let us define a general expansion 
of ff (s) as 

i S 1 // S 2 I hi . . 

+ m| + 2 A+ m| + 6 A+ m| + ---J' 

where A + , A+ and A™ are the slope, curvature and cubic expansion parameter respectively. 
Numerically, when employing the following isospin averages for the meson masses M w = 
138.0 MeV, M K = 495.7 MeV and the mass of the charged K* meson, M K * = 891.66 MeV 
[37], which is relevant for the leptonic decays of the K meson, these are found to be: 

A+ = 25.6 -HT 3 , A" = 1.31 • 10" 3 , A+ = 9.74 -lO" 5 . (12) 

The parametric uncertainties on these results are rather small. However, it is difficult to 
estimate the systematic uncertainties and therefore we have chosen not to attach errors 
to the results of eq. (12). Recent experimental results on the slope A^_ and the curvature 
A^_ have been collected in table 1. One observes that our results of eq. (12) are in nice 
agreement with the very recent measurement by KLOE [38], and in reasonable agreement 
with ISTRA [39] as well as NA48 [40], however about 2.5 a away from the KTEV result [41]. 



F™(s) = Ff "(0) 
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Collaboration A+ [1(T 3 ] A+ [10 



-3 



ISTRA 04 [39] 23.2 ± 1.6 0.84 ± 0.41 

KTEV 04 [41] 20.64 ± 1.75 3.20 ± 0.69 

NA48 04 [40] 28.0 ± 2.4 0.2 ± 0.5 

KLOE 06 [38] 25.5 ± 1.8 1.4 ±0.8 



Table 1: Recent experimental results on the slope A + and the curvature A+ of the vector 



form factor F+ n (s) in units of 10 3 . 



Since the r lepton can also decay hadronically into the second vector resonance K*' = 
if* (1410), this particle should be included in our parametrisation of the vector form factor 
F+ n (s). A parametrisation which is motivated by the RxPT framework [21,22] can be 
written as follows: 



Re[H Kn (s)+H Krl (s)\ 



m k* - s - iM K *T K ,(s) M 2 Kt , - s - iM K ,iY K ,i (s) 

(13) 



K* 

The relation of the parameter 7 to the R^PT couplings takes the form 7 = F V G V /F% — 1, 
when one assumes a vanishing form factor at large s in the Nc to infinity limit. It is 
difficult, to asses a precise value for 7, but from the work of ref. [25], we infer that it 
should be small and positive (0 < 7 1). We shall come back to the parameter 7 below. 
The width of the second resonance cannot be set unambiguously. Therefore, we have 
decided to endow the i^*(1410) contribution with a generic width as expected for a vector 
resonance. Hence, T K ^(s) will be taken to have the form 

r - (s) ^-'4:(«))- (14) 

3 The differential decay distribution 

As a last ingredient for a prediction of the differential decay distribution of the decay 
r — > v T Kii according to eq. (1), we require the scalar form factor F^ n (s). This form factor 
was calculated in a series of articles [42-44] in the framework of R^PT, again also employing 
constraints from dispersion theory as well as the short-distance behaviour, which then lead 
to a determination of the strange quark mass m s in [44]. Quite recently, the determination 
of F^ Ctc (s) was updated in [26] by employing novel experimental constraints on the form 
factor at the Callan-Treiman point Aktt, and here we shall also make use of this update. 

A remaining question is which value to use for the form factors F_F 7r (s) and F x,r (s) 
at the origin. From our chiral description of the vector form factor (5), one would obtain 
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F_^ 7r (0) = 0.978. An average over recent determinations from lattice QCD and effective 
field theory approaches [45-50], however, yields 

Ff 7r (0) = F X7r (0) = 0.972 ±0.012, (15) 

somewhat lower and also compatible with the original estimate by Leutwyler and Roos [51]. 
Furthermore, this is the value to which the scalar form factor had been normalised in 
ref. [26]. Nevertheless, what is required in the hadronic r decays is only the product 
|Kis|-F 1 +" 7r (0) which experimentally is known more precisely. Thus, in what follows, to 
normalise the form factors, we shall employ the most recent average for this number [52]: 

\V us \F^(0) = 0.2173 ±0.0008. (16) 
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Figure 2: Main result for the differential decay distribution of the decay r — > v t Ktx, 
together with the individual contributions from the K* and if* (1410) vector mesons as 
well as the scalar component residing in the scalar form factor F^ n (s). 

Our main result for the differential decay distribution is displayed as the solid line in 
figure 2, together with the individual contributions. Let us discuss our inputs and the 
individual contributions in more detail. Since the K* meson in r decays is the charged 
one, in contrast to the fit result of eq. (9), for the K* mass we have employed the PDG 
value Mjf.- = 891.66 ± 0.26 MeV [37]. However, for Gy and the barrier factor parameter 
r, the fit results of (9) have been used. The resulting contribution of the K* meson to the 
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spectral distribution is shown as the long-dashed line in figure 2. Integrating this part over 
the phase space and varying the input parameters, one finds 

B[t -> u T K*(892)} = (1.253±g;g?§ ±0.019)% = (1.253 ± 0.078)% . (17) 

The first uncertainty represents an estimate of higher order chiral corrections. To this end, 
in the exponential of eq. (5), we have replaced the factor 1/F% by l/(FfcF n ) or by 1/Fq 
with F = 87MeV being the pion decay constant at the leading order, which should give an 
idea about unaccounted further chiral corrections. The remaining uncertainty arises from 
a variation of the fit parameters Gy and r of eq. (9) and the value (16) for | V US \F^ 1T (Q). 

Next, the contribution from the scalar form factor F^ w (s) in figure 2 is displayed 
as the dotted line. Its most important contribution arises in the region below the K* 
resonance where the low-lying scalar _K"q(800) resonance is active. Integrating over the 
scalar contribution, we obtain 

B[ T - ^(X^s-wave] = (3.88 ± 0. 19) • 10~ 4 , (18) 

where the error dominantly is due to a variation of the form factor shape as discussed in 
refs. [26,43,45]. Since the scalar resonances are not well described by Breit-Wigners and 
there is also a strong interference between the dynamically generated i^o(800) and the 
pre-existing (at N c — > oo) fTg(1430) resonance, we prefer not to resolve the Kn S-wave 
contribution into individual components. (For some remarks on the _K"q(800), also known 
as the k, see section 7 of ref. [42].) The sum of the scalar and K* contributions in figure 2 
is shown as the dashed-dotted line. 

The last remaining contribution is the one due to the second vector resonance K*(1A10). 
For its mass and width, we have employed the PDG values Mk* (1410) — 1414 MeV and 
r_^*(i4io) = 232 MeV [37]. The f^*(1410) contribution turns out to depend very sensitively 
on the mixing parameter 7 defined in eq. (13), for which an estimate can be obtained on the 
basis of the total branching fraction B[r — > v t Ktt]. Adding the results of the most recent 
compilation [6], one finds B[r — > v T Kir] = (1.33 ±0.05) %. Then adjusting 7 such that the 
experimental total branching fraction, including its uncertainty, is reproduced, we obtain 
7 = 0.013 ± 0.017, in agreement with the expectation of the last section, that 7 should 
be small and positive. The contribution of the if* (1410) resonance with the central value 
of 7 is shown as the short-dashed line in figure 2. Even though this contribution appears 
rather small, because of the interference with the leading K* resonance, its influence in the 
energy region above the K* resonance is roughly as important as the scalar component. 
This is also reflected in the corresponding total branching fraction which turns out to be 

B[t -> u T (K*(892) +K*(U10))\ - B[r -> u T K*(892)} « 3.9 • 10" 4 , (19) 

however, with rather large uncertainties, but much bigger than the if *(1410) branching 
ratio by itself whose central value reads B[r — > u T K*(lA10)] = 2.1 • 10~ 6 . We notice that 
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the central result (19) is similar to the scalar branching fraction (18) and compatible to an 
estimate presented by the ALEPH collaboration [8] which yielded B[r — > u T K*(1410)] = 
(1.5lJ;o) ' 10~ 3 , but where, however, the scalar component had been neglected. 

4 Conclusions 

Upcoming much improved results on the branching fractions and differential distributions 
of hadronic r decays from the 5-factories BABAR and BELLE, as well as in the near 
future from BESIII, necessitate an analogous improvement also on the theoretical side. A 
step in this direction is taken in the work at hand, where we have presented a description 
of the decay spectrum of the decay r — > v T Kir (1) in the framework of RxPT [21,22]. 

Our approach follows the lines of an analogous description for the pion form factor 
[23-25], which employs all present knowledge of effective theories, short-distance QCD, the 
\&rge-N c expansion as well as analyticity and unitarity, and was successful in describing 
the experimental data for the pion form factor. Our central result for the Kn vector 
form factor, also including the second vector resonance if* (1410) in this channel, has been 
presented in eq. (13). As a by-product, we have determined the slope and curvature of 

\' + = 25.6 -10" 3 , A+ = 1.31 • 10" 3 , (20) 

in very good agreement with the recent KLOE results [38]. The required scalar Kir form 
factor F ( f- 7r (s) has been employed from the recent update [26], following the previous anal- 
yses [42,43]. 

Our central result for the differential decay distribution dT Kn /d v ^s of eq. (1) is displayed 
in figure 2, where also the separate contributions originating from the K*, the if* (1410), 
and the scalar component have been shown separately. 3 As can be observed from this 
figure, the S-wave Kn contribution is most prominent in the energy region below the K* 
resonance, and dominated by the dynamically generated if q (800) resonance (also known 
as the k) . Thus, by analysing experimental data in this energy range, valuable information 
about the still much debated if q (800) resonance can be obtained. Integrating the scalar 
component over the phase space, we obtain 

B[ T - ^(K^s-wave] = (3.88 ± 0. 19) • lO" 4 , (21) 

which is already rather precise, so that it will be difficult to improve this accuracy experi- 
mentally. 

The contribution of the second vector resonance if* (1410) dominantly depends on the 
mixing parameter 7. However, this parameter can be inferred from the total r — > u t Ktt 

3 A Fortran routine to generate this distribution can be obtained on request from the authors. 
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branching fraction, resulting in the estimate 7 = 0.013 ±0.017. Even though the if* (1410) 
contribution by itself is small, through interference with the K* resonance, in the range 
around 1.4 GeV it becomes about as important as the scalar resonance in this region, the 
if,* (1430). Therefore, with a dedicated experimental analysis in the region above the K* 
peak, it should be possible to also obtain information on the if* (1410). 

In summary, in this work we have started a dedicated effort to improve the description 
of exclusive strangeness-changing hadronic r-decays. We are convinced that our description 
of the decay spectrum for the decay r — > u T Kn will proof valuable in the analysis of the 
high-statistics data on hadronic r decays, acquired at the 5-factories BABAR and BELLE. 
We plan to return to this subject with other decay channels in the future. 
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